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1. INTRODUCTION 

We describe two recent additions to SFIRE, a fire spread model coupled with WRF ||Mandel et aT7||2 009 
This model builds on the earlier CAWFE code {Clark et al-1|1996| |2004| |Coen||2005) ; see (Mandel et all 



2011 



2011 



and |http: / /www, openwfm. org/wiki/WRF-Fire_development_notes for further historical details 



and acknowledgements. The coupled model is available from |0penWFM. org| An earlier version of the code is 
included in WRF release as WRF-Fire. 

2. FUEL MOISTURE MODEL 

Fire spread rate depends strongly on the moisture contents of the fuel. In fact, the spread rate drops to zero 



when the moisture reaches the so-called extinction value | |Pyne et aL||1996| . For this reason, we have coupled 



the fire spread model with a simple fuel moisture model integrated in SFIRE and run independently at points of 
the mesh. See {Nelson Jr]|2000[ |Weise et al.||2005) for other, much more sophis ticated models. (|Coen||2 005> 
used an assumed diurnal dependence of fuel moisture on time. Empirical models {Fosberg and Deemmg||l971| 
Van Wagner and Pickett| p 985 > attempt to predict fuel moisture from meteorological conditions measured daily. 



We use a simple timelag differential equation at every point of the domain. This equations has solutions which 
approach an equilibrium exponentially, if the equilibrium does not change. In general, the solutions track a changing 
equilibrium with a delay. 

First, relative humidity is computed from the atmospheric temperature T (K), waver vapor c ontents Q (kg /kg), 
and p ressure P (Pa). The saturated water vapor pressure P W s (Pa) is approximated, following (Murphy and Koop 
[2005] Eq. (10)), as P W s = exp(54.842763-6763.22/T-4.2101ogT+0.000367T+tanh {0.0415 (T — 218.8)} (53.878- 
1331.22/T- 9.44523 logT + 0.014025T)), 123 < T < 332K. The water vapor pressure is P w = PQ/ (e + ( l-e)Q), 
where e = 0.622 is the ratio of the molecular weight of water (18.02 g/mol) to molecular weight of dry air (28.97 
g/mol). We then obtain the relative humidity H = l00P w /P ws .Jhe temperature and the relative humidity of the air 
are then used to estimate the drying and wetting fuel equilibrium moisture contents {Van Wagner and Pickett|1 985 
eq. (4), (5)) {Viney|1991| eq. (7), (8)), 

E d = 0.924iT°- 679 + 0.000499e o m + 0.18(21.1 + 273.15 - T)(l - e" 0115ff ), 
E w = 0.6187J 753 + 0.000454e o m + 0.18(21.1 + 273.15 - T)(l - e -°- 115H ). 
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Figure 1 : Response of fine fuels to rain over 24 hours (a) following (Van Wa gner and Pickett| [l985) (b) from the 
time-lag model {2} by a calibration of coefficients. 

The fuel is considered as a combination of time-lag classes, and the fuel moisture contents m k (t) in each class 
k = 1, . . . , N is then modeled by the standard time-lag equation 

I ^^\1m k >E d 

1 Q\\E d <m k <E w (1) 



dt 



E ™~ k mk if m k < E u 



where T k is the lag time. We use the standard model with the fuel consisting of components with T k = 1, 10, and 100 
hour lag time, with the proportions w k > 0, Y, k =i w k = 1. given by the fuel category description i jScott and Burgan 
20051. The overall fuel moisture then is the weighted average m = J2 k =i w k m k 



)uring rain, the equilibrium moisture E d or E w is replaced by the saturation moisture contents S, and equation 
{1} is modified to achieve the rain-wetting lag time T rk for heavy rain only asymptotically, when the rain intensity r 
(mm/h) is large: 

dm k S~m k ( ( r-r \\ 

1 - exp , if r > r , (2) 



dt T rk \ \ r k 

where r is the threshold rain intensity below which no perceptible wetting occurs, and r k is the saturation rain 
intensity for fuel component k. At the saturation rain intensity, 1 — 1/e « 63% of the maximal rain-wetting rate is 
achieved. The coefficients can be calibrated to achieve a similar behavior as accepted empirical models (Fosberg 
and D eeming] |1971||Van Wagner and Pickett||1 985| >. This model estimates the fuel moisture as a function of the 



initial moisture contents and rain accumulation over 24 hours. Assuming steady rain over the 24 hours, we have 
obtained a reasonable match with the Canadian fire danger rating system (Van Wagn er and Pickett||1985) using 
S = 250%, T rk = 14 h, r Q = 0.05 mm/h and r k — 8 mm/h, cf., Fig. QJ Of course, since fuel may dry up 
between episodes of rain, in the case of intermittent rain, the result of the model depends also on the temporal 
rain distribution, not only the total accumulation - just like fuel moisture in reality. 

Because we want the model to support an arbitrarily long time step, an adaptive exponential method was 
implemented. The method is exact for long time step when the atmospheric variables and the rain intensity are 
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Figure 2: (a) Simulation of fuel moisture contents at a single point for fine fuels. The drying and wetting equilibria are 
computed from the WRF state. The fuel moisture contents does not change when it is between the two equilibria. 
The red vertical lines correspond to periods of rain, where the equilibrium is 2.5 (above the range shown). The fuel 
moisture contents increases during rain and exhibits diurnal variability. The flat part between 100 and 200 hours is 
due to one day of missing data, (b) Advection of smoke by tracers injected in WRF-Chem, in an ideal example. 



constant in time. Equations {T} and |2) have the common form 

dm E — m 
~dt = T ' 



(3) 



For longer time steps, use the exact solution of {3} over the time step interfal [t n ,t n+ i], with constant coefficients 
taken as their values at t n+1/2 = t n + At/2, which gives M n+1 = M n + (E n+1/2 - M n ) (l - e -A */ T »+v 3 ), 



where E. 



n+l/2 



and T n 



-1/2 



are computed at t, 



-1/2- 



This method is of second order and it is particularly useful 



when the time step is comparable to or even larger than the time lag and the coefficients E and T vary slowly, 
such as in the case of fast drying just after the rain ends. For very short time steps, however, the rounding 
error in the subtraction of the almost equal quantities in 1 - e - At / T «+i/2 will pollute the solution. Thus, for 
At/T n+1/2 < e — 0.01, we replace the exponential by a truncated Taylor expansion, and use the second order 



method M n+1 = M n + (E. 



ra+l/2 



M n ) 



At 



1 



A/ 



Because the time scale of the moisture changes (hours) is very different from the time scales of the atmoshere 
(minutes( and fire (seconds), the moisture model runs at a multiple of the WRF time step. WRF variables p 
(pressure, Pa), T2 (temperature at 2m, K), and Q2 (water vapor mixing ratio, kg/kg) at the beginning and the end of 
the moisture model time step are averaged to obtain the values of P, T, and Q used in the timestep of the moisture 
model, 

P(t n+ i)+P(i„) _ T2 (t n+ i) + T2 (t n ) _Q2(t n+1 ) + Q2(t n ) 

— 7) -i J-n+1/2 — 7) j Wnn+1/2 — 7y > 
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This is done for 2nd order accuracy as well as for compatibility with the computation of the rain intensity from the 
difference of the accumulated rain, 

_ RAINC(i„+i) + RAINNC (i„+i) - RAINC(t„) - RAINNC (t n ) 
r n+ i/2 - ^ • 

These values are then used to compute the equilibrium moisture contents and the time lag, and one time step is 
performed, as described above. See Fig. [2](a) for a typical simulation result. 

The fire model runs on a finer grid than the atmospheric model, typically refined by a factor of 10 or more. The 
moisture model consists of two steps. 

In the first step, we compute the moisture content m k of the fuel components on the nodes of the atmospheric 
grid on the Earth surface for several reasons. (1) The WRF variables are known at the nodes of the atmospheric 
grid and no interpolation is needed. (2) The atmospheric grid is relatively coarse, so the added costs of storage 
of several surface moisture fields m k and of the computation are not significant. (3) The computation is done for 
each fuel component indivudually over the whole domain and it does not depend on the actual fuel map. (4) The 
fire model does not need to run for the computation of m k . 

In the second step, we interpolate the values of m k to the fire grid and compute the weighted averages on the 
fire grid following using the actual fuel map. There are several options how to run the fuel moisture model: (1) 
Turn on both steps of the moisture model: compute the moisture fields m k and interpolate on the fire grid, as the 
fire model runs. This option is intended for the actual fire simulation. (2) Turn on the first step, computation of the 
moisture fields m k only, for an extended run (many days) to evolve the m k in response to the simulated weather, into 
a dynamic equilibrium. The moisture fields are stored in WRF state in the output and the restart files. (3) Run the 
first step, computation of the moisture fields m k , as a standalone executable from stored output files from standard 
WRF runs, adding the moisture fields m k to the files. (4) Start the coupled atmosphere-fire-moisture simulation 
from the moisture fields evolved over time as above. (5) Run the fire model as a standalone exectutable, without 
feedback on the atmosphere, with fuel moisture computed using the second step only from the fields m k in WRF 
output files produced in advance. 

3. COUPLING WITH WRF-CHEM 

Coupling with WRF-Chem is implemented by inserting the smoke intensity in the WRF-Chem arrays emis.ant 
and tracer at the ground layer and with the species index p.smoke. WRF needs to be configured with the chem 
option and built with em_f ire or em_reai, and the appropriate options trace_opt=i and chem_opt=i4 set in the 
namelist. A sample visualization of an ideal run with smoke transport is in Fig.[2jb). See the WRF SFIRE Users' 
Guide at |http : //www, openwfm. org/wiki/Users_guide| for further details and more information as the code 
develops further. 

Currently, the smoke inserted into WRF-Chem is simply proportional to the fire heat flux. However, a more 
sophisticated scheme following the FEPS (Fire Emission Production Simulator) is being implemented (And erson 
|et al.||2004) . In this new method, the fire emission will be treated as the fluxes of carbon dioxide (C02), carbon 
monoxide (CO), methane (CH4) and particulate matter (PM 2.5) released into the atmosphere at the location of 
the fire. The total emission for each of these species will be computed based on the fuel burn rate (amount of fuel 
burnt per unit time), fuel load, and the smoldering correction computed based from the wind speed and relative 
humidity. The PM 2.5 will be portioned into the the accumulation and nuclei mode and released into the atmosphere 
as PM25J and PM25I respectively. 

Aside from the obvious application to smoke dispersal, which is important in practice, taking into account the 
composition of the fire emitted smoke will allow for treating it as chemically and physically active. The species 
released into the atmosphere by the fire will undergo chemical reactions allowing for capturing the secondary aerosol 
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effects. The particular mater will interact with atmospheric radiation, and have a potential to serve as condensation 
nuclei when suitable cloud microphysics is available. This may be of great importance in the simulation of visible 
plumes and pyrocumulus clouds, under development in | |Peace et al.|201 1) . 
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